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The effects of the number of significant figures (NSF) in the interpolation polynomial coefficients (IPCs) of the 


weighted sum of gray gases model (WSGM) on results of numerical investigations and WSGM optimization were 


investigated. The investigation was conducted using numerical simulations of the processes inside a pulverized 


coal-fired furnace. The radiative properties of the gas phase were determined using the simple gray gas model 
(SG), two-term WSGM (W2), and three-term WSGM (W3). Ten sets of the IPCs with the same NSF were formed 
for every weighting coefficient in both W2 and W3. The average and maximal relative difference values of the 


flame temperatures, wall temperatures, and wall heat fluxes were determined. The investigation showed that the 


results of numerical investigations were affected by the NSF unless it exceeded certain value. The increase in the 
NSF did not necessarily lead to WSGM optimization. The combination of the NSF (CNSF) was the necessary 


requirement for WSGM optimization. 
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Introduction 


In numerical investigations of the processes taking 
place inside pulverized coal-fired furnaces, the radiative 
properties of the gas-phase combustion products are 
often modeled using the weighted sum of gray gases 
model (WSGM) [1, 2], by which the real gas is replaced 
by a mixture consisting of several gray gases plus one 
clear gas. The radiative properties of every gray gas are 
defined by a temperature-dependent weighting coeffi- 
cient and constant absorption coefficient. The form of the 
weighting coefficient dependence on temperature is a 
polynomial of selected order. 

Several WSGMs have been developed so far, inclu- 
ding spectral line based WSGM [3], narrow band based 
WSGM [4], and WSGM based on total radiative pro- 
perties [5]. The effects of the number of gray gases in the 


Nomenclature 

A surface area [m’] 

a weighting coefficient [—] 
B polynomial coefficient [—] 


E, blackbody emissive power [W m7] 


WSGM on the calculation of the radiative heat transfer 
for the spectral line based WSGM was described by Coe- 
lho [6], who used 3 and 20 gray gases to evaluate perfor- 
mances of two radiation models for the calculation of 
wall fluxes and radiative heat sources inside a rectangular 
3-D enclosure with fixed temperature and concentration 
distributions. The exact solution was obtained using the 
statistical narrow-band model [7]. An increase in the 
number of gray gases from 3 to 20 reduced the average 
and maximal relative errors, but in many cases the level 
of reduction was small. Park and Kim [8] used the re- 
grouping technique to investigate the possibility of redu- 
cing the number of gray gases of the narrow band based 
WSGM. The exact values of the wall fluxes and average 
intensity of radiation inside the cubic 3-D enclosure with 
fixed temperature and concentration spatial distributions 
were found using the statistical narrow-band model. 


dissipation of the turbulent kinetic energy 


E [m? s] 

p density [kg m”] 

@P time-averaged general variable 
Subscripts 


G;G; volume-volume total exchange area [m7] 
GG; volume—volume directed flux area [m] 
GS; volume-surface total exchange area [m] 
GS; volume-surface directed flux area [m’] 
H enthalpy [J kg"'] 

H convection heat transfer coefficient [W m? K] 
k turbulent kinetic energy [m s] 

k; thermal conductivity [W m" K"] 

K, absorption coefficient [m"'] 

K; scattering coefficient [m'] 

l layer thickness [m] 

N total number [—] 

q heat flux [W m7] 

Ò heat transfer rate [W] 

S/S; surface—surface total exchange area [m°] 
SS; surface—surface directed flux area [m?] 
SiG; surface—volume directed flux area [m] 
So source term for the gas-phase variable Ø 
T temperature [K] 

U time-averaged velocity vector [ms”'] 

V volume [m°] 

Greek symbols 

T transport coefficient [kg m” s] 

ô relative difference [%] 

ò total emissivity [-] 


Nine hundred gray gases were replaced by 3, 5, 7, 10, 15, 
and 20 gray gases. The average and maximal relative 
errors generally decreased with the increase in the 
number of gray gases and when the number of gray gases 
became bigger than 5, the reduction of the average and 
maximal relative errors became small. Investigations 
reported in [7, 8] showed that the increase in the number 
of gray gases in WSGM led to exact solution, and also 
showed that the average and maximal errors were not 
considerably reduced when the number of gray gases 
exceeded a certain number. That results, obtained for 
enclosures with fixed temperature and concentration 
fields, indicated the possibilities of the optimization of 
the number of gray gases in WSGM, that is WSGM 
optimization. 

Crnomarkovic et al. [9] investigated the optimization 
of WSGM based on total radiative properties. 
Investigations were carried out for the process inside a 
pulverized coal-fired furnace of a utility scale boiler. All 


al ash layer 

bw boiling water 

dp dispersed phase 

gp gas phase 

88 gray gas 

f flame 

ml metal layer 

p particle 

r total number of significant figures 

rad radiation 

SZ surface zones 

VZ volume zones 

w wall 

D related to variable @ 

Abbreviations 

CNSF combination of the number of significant 
figures 

IPC interpolation polynomial coefficient 

NSF number of significant figures 

RPM radiative property model 

SG simple gray gas model 

TEA total exchange area 

WSGM weighted sum of gray gases model 

W2 two-terms WSGM 

W3 three-terms WSGM 


variables of the flow field were obtained from the 
numerical simulations. Gas-phase radiative properties 
were determined by the simple gray gas model (SG), 
two-term (one gray gas plus one clear gas) WSGM (W2), 
and three-term (two gray gases plus one clear gas) 
WSGM (W3). Three sets of the flame radiative properties 
were formed. The interpolation polynomial coefficients 
(IPCs) were determined with seven significant figures. 
The investigation showed that WSGM could be opti- 
mized for one set of the flame radiative properties. 

This research was conducted to further elucidate the 
effects of the number of significant figures (NSF) in the 
IPCs on numerical investigations. Objectives of the in- 
vestigation were twofold: 

(1) To find effects of the NSF on results of numerical 
investigations. This research was expected to reveal 
possible limitation in the NSF. 

(2) To find effects of the NSF on WSGM optimization. 
This research was expected to show if the increase in the 
NSF necessarily led to WSGM optimization. 


This research used the gas-phase radiative property 
models from previous one [9]. Ten sets of the IPCs with 
the same NSF for each weighting coefficient in W2 and 
W3 were formed and numerical simulations were run for 
every set of the IPC. As the exact values of variables 
were not known, the relative differences of the selected 
variables were used instead of relative errors. The term 
WSGM optimization implied the decrease of the relative 
differences of variables with the increase in the number 
of gray gases. 

The investigation was conducted using an in-house 
developed computer code. The shape and dimensions of 
the case-study furnace and coal properties are described 
in [10]. The furnace of a 210 MW monoblock unit is 
tangentially fired with domestic lignite Kolubara. The 
furnace is 45.0 m high, 16.5 m wide, and 14.5 m deep. 
The furnace is equipped with six jet burners, five of 
which being in constant operation in the test-cases con- 
sidered. The following three sections outline the mathe- 
matical model with the method of analysis, obtained 
results with discussion, and conclusions. 


Mathematical model and determination of the 
relative differences 


A description of the mathematical model of the reac- 
tive turbulent two-phase flow with radiative heat ex- 
change process taking place inside the furnace can be 
found in [9-12]. Here, the most important model chara- 
cteristics are described. 

The gas phase is described by the time-averaged 
partial differential equations for the conservations of 
momentum, enthalpy, gas-phase concentrations, turbulent 
kinetic energy, and its dissipation rate in the Eulerian 
reference frame. The general form of the gas-phase 
equation is the following: 


div(pU®) = div(Pggrad®) +S +So, (1) 


where Sop is the source term for ® due to presence of 
particles. The set of equations is closed by k-e turbulence 
model. The pressure field is solved using the SIMPLE 
algorithm. The flame temperature distribution is obtained 
for the condition of thermodynamic equilibrium between 
the gas and dispersed phases [11]. 

The dispersed phase is described by the differential 
equations for motion and changes in mass and energy in 
the Lagrangian reference frame, with PSI-Cell method 
for influence on the gas phase. Heterogeneous reactions 
taking place during coal combustion are modeled in the 
kinetic-diffusion regime, with the combustion model as 
described in [12]. 

Radiative heat exchange is determined by the Hottel’s 
zonal model, by which furnace space and walls are 


divided into volume and surface zones, respectively. The 
purpose of radiation model in numerical simulations is to 
find the source term due to radiation for the enthalpy 
equation as well as wall heat fluxes. When WSGM is 
used, the source term for the enthalpy equation is the 


following [13]: 
Nyz Noz 
SH radi = > GG Eom + DS GEbn 
n n=l (2) 
-49 ag (T)K a ViEp js i=1,...,N,, 
j=0 
where 


The index j = 0 designates the clear gas. TEAs are 
determined for the radiative properties of the gray gases. 

The wall heat flux of a surface zone s; is the difference 
between energy absorbed and lost divided by the surface 
area. When WSGM is used, it is expressed as: 


N, N 
> GSEs + DSS Eon E AOE, ; 
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Heat transfer rate through a surface zone s; is a pro- 
duct of the wall heat flux and surface area: Q, = Iw iA; » 


whereas the rate at which heat is absorbed by boiling 
water is a summation of a heat transfer rates through all 
surface zones which represented the solid wall: 


Ow a >ü : 

When SG is used, the directed flux areas in eqs. (2) 
and (3) are replaced by TEAs determined for one set of 
the flame radiative properties (absorption and scattering 
coefficients). 

The absorption coefficients of the gas phase were ca- 
Iculated from total emissivities data and for temperature 
T = 1500 K. The gas-phase emissivities were found in 
temperature range of 600-2400 K with a step of 50 K, 
using Leckner’s model [14]. The absorption coefficient 
for the SG was Ky.) = 0.07 m”. The absorption coe- 
fficients for the W2 and W3 were determined by the 
method described in [5]. For a two-phase mixture of gas 
with particles, each gray gas in WSGM actually becomes 
a gray medium whose absorption coefficient is the sum 


of absorption coefficients of the gray gas and dispersed 
phase and whose scattering coefficient is that of the 
dispersed phase [15]. The radiative properties of the gray 
gases for the W2 and W3 are given in Table 1. 


Table 1 Radiative properties for the W2 and W3 


W2 W3 
88 Ka K, Ka Ky 
0 1.0 1.0 1.0 1.0 
1 1.23 1.0 1.014 1.0 
2 - - 1.41 1.0 


The weighting coefficients were determined by fitting 
the emissivities data [5]. On the basis of previous 
investigation [9], the values of the weighting coefficients 
were interpolated using the third order polynomial: 


3 r V 
a,(T)= ) b, ;| —— 4 
z( ) 2 me ( ) 

The coefficients b,; were determined with 10 
significant figures via the least-squares method [16]. The 
coefficients are presented in Table 2 in normalized 
decimal floating-point form: +0.dd)...d,x10” [17], 

The integer n > 0, as well as 1 <d,<9 and 0 < d} < 9 

. Neg 

for each k = 2,..,r. The coefficient ay =1.0-)) 24, . 
The subscript r designates the total NSF [18] and each 
digit d, (k = 1,...,r) is a significant figure. The NSF in 
every coefficient was reduced by one using the rounding 
method until the polynomial coefficients with the single 
NSF were obtained. 

The furnace walls consisted of two layers: a 4.0 mm 
thick metal wall layer and a 0.3 mm thick ash deposit 
layer. Wall temperatures are determined on the basis of 
qw and conductive and convective heat transfer: 


Tas =T tana pth | (5) 
f "(A k, mi kial 

where Tew = 615 K. The value of the convective heat 
transfer coefficient (h = 14.0 kW m? K`’) on the side of 
boiling water was adopted from [19]. The thermal con- 
ductivity of ash deposit layer is actually the effective 
thermal conductivity and its dependence on temperature 
is that from [20] (the sample whose silica ratio is 62), 
whereas the thermal conductivity of the metal layer is 
that of carbon steel [21, 22]. Wall emissivity is 0.8. The 
gas-phase thermophysical properties are determined from 
the equation of state, tabulated values, and empirical rela- 
ions. The set of equations is solved by the finite-diffe- 
rence method. Discretization and linearization of the 
gas-phase equations are achieved using the control 
volumes and hybrid difference scheme methods. Iterative 
procedure stability is provided by the under-relaxation 


method. 


Table 2 Coefficients in the interpolation polynomials 


W2 W3 
j 
dı dı da2 
0 0.4246901182 0.3659433153 0.5638441144x10"! 
1 0.4030451291 0.4457903327 —0.1790604993 
2  —0.2618388217 —0.3326471831 0.2368749224 
3 0.4089782078x10' 0.5355003890x107 -0.4040215113x10" 


The relative differences are determined for the fol- 
lowing variables: flame temperature 7;, wall heat flux qy, 
and wall temperature 7,,. To find effects of the NSF on 
results of numerical investigations, the relative differ- 
ences were determined taking the values of the variables 
obtained with 10 significant figures as the reference 
values: 

= Paea 90.0% (6) 
hoi 

where 77 stands for a variable. The subscripts r and 10 
in eq. (6) designate the NSF in the IPCs for which the 
variables were determined. Equation (6) was used to de- 
termine the relative differences of the rate of the ab- 
sorbed heat. 

To find effects of the NSF on WSGM optimization, 
the relative differences in the considered variables be- 
tween W3 (on one side) and SG or W2 (on the other side) 
were determined taking the values obtained by the W3 as 
the reference values: 


ne Pees "199.004 (7) 


3,1 

where the subscript RPM designates SG or W2. The 
relative differences were determined for the same and 
different NSF in the W2 and W3. The pair of the NSF 
used in both W2 and W3 composed the combination of 
the NSF (CNSF). 

Average values of the relative differences were deter- 
mined as the arithmetic mean of the values over all con- 
trol volumes (for Tà and all surface zones (for qw and Ty). 
For WSGM optimization, the investigation included both 
average and maximal values of relative differences and 
also relative differences of the rate of absorbed heat. 


Results and discussion 


Radiative heat exchange was solved on the coarse nu- 
merical grid composed of the cubic volume zones with 
edge dimensions of 1.0 m. The total numbers of volume 
and surface zones were 7956 and 2712, respectively. 
Direct exchange areas of the close zones were deter- 
mined using the correlations given in [13]. TEAs were 


determined using the original emitters of radiation 
method [6], with the values of the TEAs then improved 
via the generalized Lawson’s smoothing method [23]. 
The flow field was determined on the fine numerical grid, 
which was obtained by dividing each volume zone into 
64 control volumes. The radiative energy source term of 
enthalpy equation, determined for every volume zone, 
was divided equally to all of its 64 control volumes. The 
total number of control volumes was 620 136. Agreement 
with experimental data and grid independence study were 
proved in [11], with the gas-phase radiative properties 
determined by the SG model. All results were obtained 
after 4000 iterations. 

The effects of the NSF on the results of numerical in- 
vestigations through the average values of the relative 
differences in the variables are illustrated in Table 3. The 
results are not changed when the NSF is greater than 9 
for W2 and when the NSF is greater than 8 for W3. Thus, 
totally 9 x 8 = 72 combinations should be used to ex- 
amine WSG optimization. For that, it is necessary to find 
all CNSFs for which the average and maximal relative 
differences decrease with the increase in the number of 
gray gases. 


Table 3 Effects of the NSF on average values of relative dif- 
ferences [%] 


qw Tw Te 
i w2 W3 w2 W3 w2 w3 
1 1.37 3.07 0.21 0.51 0.43 0.37 
0.91 1.72 0.13 0.27 0.27 0.36 
0.93 1.11 0.13 0.16 0.28 0.40 
0.98 0.85 0.15 0.12 0.36 0.26 
1.11 1.55 0.15 0.20 0.36 0.46 
0.88 1.09 0.13 0.14 0.25 0.35 


0.72 0.86 0.10 0.12 0.20 0.23 
1.02 0.00 0.14 0.00 0.39 0.00 
0.00 0.00 0.00 0.00 0.00 0.00 


Oo oo u Dn A U N 


Among 72 combinations, three of them satisfy the re- 
quirement for WSGM optimization, and they are shown 
in Tables 4 and 5. Table 6 shows relative differences of 
the rate of the absorbed heat. 

Every CNSF is represented by a pair r,r, where the 
first r designates the NSF in the W2 and the second r 
designates the NSF in the W3. The results presented in 
Tables 4—6 show that the best results is obtained for 
CNSF 7,7. The average value of the relative difference 


Osc indicates the agreement of the results obtained for 
W3 with 7 significant figures with the experimental data. 
The average value of the relative difference dw2 shows 
that the values of the selected variables are changed very 
little with the increase in the number of gray gases. On 
the other hand, calculation time is shorter by about 24 
hours if W2 is used instead of W3 [9]. On the basis of 
this investigation, application of the W2 model with the 
polynomial coefficients determined with 7 significant 
figures for numerical investigations would be optimal. 

The increase in the NSF in both W2 and W3 does not 
necessarily lead to WSGM optimization. CNSF for 
which the number of gray gases in WSGM can be opti- 
mized must be found by numerical simulations and 
comparing the relative differences. Until establishing 
some shorter procedure, the only way to find the CNSF 
for WSGM optimization is as outlined in this research. 

Figs. 1-3 show the distribution of the selected va- 
riables and their relative differences for CNSF 7,7. Fig. 
la—c shows the flame temperatures in the vertical middle 
cross-section and the corresponding relative differences. 
The distribution of the relative difference ðw resembles 
that of the dsg with the reduced values of the relative 
differences. Figs. 2a—c and 3a—c show the wall fluxes and 
surface temperatures over the rear furnace wall with the 
corresponding relative differences. The relative differ- 
ences dw. are much smaller than ósa so that no resem- 
blance is found. 

This investigation shows that the NSF in the poly- 
nomial coefficients affects the results of numerical inves- 
tigations. For the selected furnace and processes inside it, 
WSGMs with more than three terms (two gray gases plus 
one clear gas) are not needed. After the optimization is 
proved, it is enough to use two-terms WSGM. That 
number of gray gases is smaller than proposed for radia- 
tion of the gaseous medium [24, 25], probably because of 
the presence of dispersed gray particles, which contribute 
to the grayness of the medium. It does not mean that 
WSGM with more than three terms is not needed for 
other pulverized coal-fired furnaces under different oper- 
ating conditions. Whenever the WSGM is used, the nu- 
mber of gray gases should be somehow justified. If WSGM 
optimization cannot be proved, number of gray gases 
should be selected using other criteria or assumptions. In 
that case, the number of gray gases could be chosen on 
the basis of the agreement with experimental data. 
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Fig.1 The flame temperature for CNSF 7,7 in the vertical middle cross-section, y = const.: a) Tẹ obtained for the W3, b) dsg, c) dw. 


a b 


Fig.2 The wall flux for CNSF 7,7 over the rear furnace wall: a) qw obtained for the W3, b) dgg, c) Owe. 


a b 


c 


Fig.3 The wall temperature for CNSF 7,7 over the rear furnace wall: a) Tw obtained for the W3, b) dg, c) ôw2- 


The zonal model of thermal radiation with homoge- 
neous radiative properties of the flame is rarely used in 
numerical simulations. Instead of it, discrete ordinate 
model, model of spherical harmonics and discrete trans- 


fer model are used. All these models were developed 
with some assumptions, which simplify the method of 
Table 4 The average values of the relative differences [%] for 
the selected CNSFs 


nr 
SG W2 SG W2 SG W2 


7,4 2.46 0.82 0.38 0.11 0.44 0.23 
TT 2.19 0.61 0.36 0.08 0.36 0.19 
9,7 2.19 0.60 0.36 0.08 0.36 0.19 


Table 5 The maximal values of the relative differences [%] 
for the selected CNSFs 


qw Tw Te 
ur 
SG W2 SG W2 SG W2 
7,4 12.24 9.94 1.93 1.24 56.47 48.85 
Tal 9.36 6.15 2.16 1.39 56.49 25.54 


9,7 9.36 4.15 2.16 0.63 56.49 25.49 


Table 6 Relative differences [%] of the rate of absorbed heat 
for the selected CNSFs 


rr OG Ow2 

7,4 1.59 0.04 
BT 1.63 0.007 
9,7 1.63 0.05 


radiative transfer solution. This research was conducted 
with zonal model to exclude the influences of the model 
assumptions on the results. Nevertheless, the procedure 
described in this paper can be used to prove WSGM op- 
timization with some other radiation model and with 
nonhomogeneous radiative properties. 


Conclusions 


In order to find the effects of the NSF in IPCs of the 
weighting coefficients on results of numerical investiga- 
tions and WSGM optimization, the computer code that 
solved the set of equations of the mathematical model 
describing the process inside the pulverized coal-fired 
furnace was used. Three gas-phase radiative properties 
models were employed: SG, W2, and W3. For W2 and 


W3, ten sets of the IPCs with the same NSF were formed. 


The numerical simulations were run for every set of IPCs 
and the relative differences in the selected variables were 
found. The conclusions can be generalized in the follow- 
ing way: 

- When the gas-phase radiative properties are deter- 
mined by WSGM, there is a specific NSF over which the 
results of numerical investigations of processes inside 
pulverized coal-fired furnaces are not changed further. 
Values of the average relative difference do not gradually 
decrease with an increase in NSF. 

- CNSF is the necessary requirement for WSGM op- 
timization. The increase in the NSF does not necessarily 
lead to WSGM optimization. CNSFs for which the num- 


ber of gray gases in WSGM can be optimized must be 
found by numerical simulations in every individual case, 
which is a time-consuming procedure. Further work is 
needed to find more efficient numerical procedure for the 
task. 
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